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0^ Abstract. We use type la supernovae (SN la) data in combination with recent baryonic 

Q acoustic oscillations (BAO) and cosmic microwave background (CMB) observations to con- 

jH strain a kink-like parametrization of the deceleration parameter (q) . This g-parametrization 

^ can be written in terms of the initial (qi) and present (go) values of the deceleration pa- 

I— I rameter, the redshift of the cosmic transition from deceleration to acceleration (zt) and the 

I redshift width of such transition (r). By assuming a flat space geometry, qi = 1/2 and adopt- 

^ ing a likelihood approach to deal with the SN la data we obtain, at the 68% confidence level 

<^ (C.L.), that: zt = 0.56tojo, ^ = ^-^'^-ofo and qo = -0.31+[!:}J when we combine BAO/CMB 

fvj observations with SN la data processed with the MLCS2k2 light-curve fitter. When in 

^ ^ - - - ' 0.13 

0.07' 



C^ this combination we use the SALT2 fitter we get instead, at the same C.L.: zt = 0.64"'' ' 



CO T = 0.361q'^7 and qo = — 0.53lg'^3. Our results indicate, with a quite general and model inde- 

^^ pendent approach, that MLCS2k2 favors Dvali-Gabadadze-Porrati-like cosmological models, 

^_^ while SALT2 favors ACDM-like ones. Progress in determining the transition redshift and/or 

L| the present value of the deceleration parameter depends crucially on solving the issue of the 

•^H difference obtained when using these two light-curve fitters. 

X 
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1 Introduction 

By the end of the last century, type la supernovae (SN la) observations [1] suggested that 
the expansion of the Universe is not slowing down, as would be expected from attractive 
gravity, but is speeding up. During the last years, in spite of several sources of systematics, 
that are still not completely well understood [2], cosmic acceleration has been confirmed 
by subsequent SN la surveys (for recent compilations see [3-6]) and also by other probes, 
like, for instance, cosmic microwave background (CMB) and baryon acoustic oscillations 
(BAO) [7, 8], that in combination have corroborated the SN la result. Discovering the source 
of cosmic acceleration is one of the biggest challenges of modern cosmology. Fortunately, 
several experiments are underway, and others are being planned, that could bring a major 
progress in the determination of the cosmic expansion history, helping us to understand the 
fundamental physics behind cosmic acceleration. 

Many theoretical approaches to investigate cosmic acceleration have been suggested. 
Here we are interested in the so-called kinematical approach [9, 10] in which the deceleration 
parameter (q) is parametrized as a function of the redshift (z). Why q{z) is considered and 
not another quantity like the equation of state parameter {w)7 The main reason for this 
choice is that with the g-parametrization only a metric theory of gravity is assumed (one is 
not restricted to general relativity) and the assumptions on the dark sector are minimum. Of 



course we should keep in mind that by parametrizing q (or any other cosmological quantity) 
one is not driven solely by the data. In fact, by defining how the deceleration parameter 
changes with redshift, a theoretical choice is performed. This choice should be considered a 
good one, if it is able to generalize, or at least describe, a large number of viable cosmological 
models at as large as possible a redshift range. By viable cosmology we mean one that 
has a kind of early matter dominated phase, in which structure might form, and a recent 
acceleration phase in accordance with observations. 

In this work we use recent SN la, BAO and CMB data to investigate the transition 
from cosmic deceleration to acceleration. We are particularly interested in the redshift of 
this transition, its (redshift) duration and the present value of the deceleration parameter. 
We remark that although the numerical values we shall obtain are conceivably sensitive to 
our specific q{z) parameterization, differently from other g-parametrizations in the literature, 
our kink-like q{z) expression is quite general and valid for a wide redshift range. As will be 
discussed in the next section, it generalizes several cosmological models and can be extended 
to high redshift. 

This paper is organized as follows: in section 2, we present our kink-like g-parametrization, 
discuss some of its properties and show how it generalizes known cosmological models. In 
section 3, we present the cosmological tests (SN la and BAO/CMB) we use to constrain the 
model parameters. In section 4 our main results are exhibited and we conclude in section 5. 

2 The kink-like q{z) parametrization 

2.1 Basic equations 

In this work, we investigate flat cosmological models whose deceleration parameter, after 
radiation domination, can be described by the following expression [10] 

9(.)=* + ^?V^. (2.1) 

If \ 1+^ ; 

Here the parameter zt denotes the transition redshift {q{zt) = 0) from cosmic deceleration 
{q > 0) to acceleration (g < 0), gj > and qj < are, respectively, the initial (z » zt) 
and final (z — t- —1) asymptotic values of the deceleration parameter. The parameter r is 
associated with the width of the transition. It is related to the derivative of q with respect 
to the redshift at z = zt, that is to the jerk (j) (see eq.(2.9) for definition) at the transition 
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In figure 1 (upper panel) we plot q{z) for zj = 1 and different values of the parameter r. In 
the lower panel we show q{z) for r = 1/3 and different values of the parameter zt- In both 
panels it has been assumed qi = 1/2 and qj = —1. 
From the definitions of the Hubble parameter. 
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Figure 1. Dependence of the g-parametrization on the parameters. Toy. dependence on r for zt — \. 
Bottom: dependence on zt for r — 1/3. In both cases qi — 0.5 and g/ = —1. 



we can write 



H = Hoexp {l + q{x))dln{l + x). 

Jx=0 



(2.5) 



In the above equations, as usual, a = 1/(1 + 2;) is the scale factor of the Friedman- Robertson- 
Walker metric, Hq = lOO/i km/(s Mpc) is the present value of the Hubble parameter and 
a dot over a quantity denotes differentiation with respect to the cosmic time. By using eq. 
(2.1), eq. (2.5) can be integrated to give 
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where, in the last equality, we have introduced the quantity 
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The above definition is very convenient and facilitates the comparison of special cases of our 
g-parametrization with more traditional cosmological models like ACDM, u^CDM, etc. At 




Figure 2. Dependence of the jerk j{z) on t for zt — 1, qi = 0.5 and qj = —1. 

this point we stress that the purpose of our g-parametrization, eq. (2.1), is to describe the 
late time transition from a cosmic decelerated to an accelerated phase. Although its validity 
can be extended to high redshift (z < 1 — 2 x 10^), it does not describe the behavior of the 
deceleration parameter in the very early universe, that is, during the radiation dominated 
era (RDE) (when q = 1). Therefore, in eq. (2.8) the limit should be understood as z >> zt, 
but after RDE. 

Another useful dimensionless quantity is the jerk [11] 






(2.9) 



With the g-parametrization given by eq. (2.1), we obtain for the jerk the following expression 
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Notice that j{zt) 
qi = 1/2, qf = -1 
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in accordance with eq. (2.2). 



(2.10) 
In figure 2 we plot j{z) for 



Zt = 1 and different values of r. Models with r > 1/3 have j{z) bellow 
one, while those with r < 1/3, have j{z) above the line j = 1. Changing zt does not alter 
this feature; the only difference is the redshift position of the maximum (r < 1/3) or the 
minimum (r > 1/3). Higher zt, higher the redshift of the j{z) maximum/minimum. 

From eq. (2.1) we can relate the final value of the deceleration parameter, qf := q(z = 



-1), to its current value qo := q{z = 0) 
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By substituting eq. (2.11) in (2.1) we get 

q{z) = — > ( ^ . (2.12) 

As observed above, with our parametrization we aim to describe the transition from a cosmic 
decelerated phase to an accelerated one. Therefore, since qi > and qf < 0, by using eq. 
(2.11) it is straightforward to show that the parameter r is constrained to the interval 

ln[l + zt] 
^<^<^rrfY (2.13) 

2.2 Special Cases 

In most scenarios, large scale structure formation requires that at early times the universe 
passes through a kind of matter dominated era in which, H^ oc (1 + z)^ and, consequently, 
q = 1/2. As discussed above, our g-parametrization is designed to describe cosmic evolution 
starting from a decelerated phase on and, for the sake of simplicity, we fix qi = 1/2 in this 
work. In fact, with this assumption we do not lose much generality, and are essentially 
excluding from our description models with a coupling in the dark sector (see reference [10] 
for more on this point). 

We present below some special cases that can be described by our parametrization. 

1. The flat standard cosmological model with constant dark energy equation of state 
(wCDM) (w := p/p = constant, where p is dark energy pressure and p its energy 
density) is recovered if we identify fimoo = ^mo and —3w = 1/r = (1 — 2qf). For 
wCDM the transition redshift is given by zt = [{0,mo — 1)(1 + 3tf)/i7mo]~ "" — 1- 
In particular, for flat cosmological constant + Cold Dark Matter cosmological model 
(ACDM) {w = -1), qf = -1, r = 1/3 and zt = [2(1 - Qruo)/ ^mo]^^^ - 1- Note also 
that for ACDM j = 1. 

2. If we apply the definition of r, given by eq. (2.2), with q, = 1/2 and qf = —1, to 
the flat DGP (DGP) brane- world model [12] we obtain r = 1/2, independent of ^mo- 
Therefore g-models with zt = (2(1 — flmo)'^/ ^mo)"^ ~ 1 (the DGP redshift transition) 
and r ~ 1/2 are expected to be a good approximation for DGP models. 

3. Neglecting baryons the flat quartessence Chaplygin model {p = — M^^^+^Yyo") (see 
[13], [14], [15] and [16]), is obtained if (besides qi = 1/2) we assume qf = —1, identify 
1/r = 3(1 + a) and Qmoo = (1 - woy/'^^+°'\ where wq = -(M^/po)"'^^ is the present 
value of the equation of state parameter. For this model zt = {—2wq/ (1 + u^o))^- 

4. Modified Polytropic Cardassian model [17], [18]. This model depends on three param- 
eters: m (denoted by q in [17]), n and Vtmo- The kinematics of this model is obtained 
if we identify, Vimo = ^moo, fn = 1/(t(1 — 2gj)) and n = 2/3(1 + qf). 



3 Observational Data 

In the next sections we derive the constraints set by three different experiments upon the 
parameters zt, r and qj (or qq). We use the following datasets: a) observations of type 
la supernovae from the first-year SDSS-II Supernova Survey together with other supernova 
datasets (ESSENCE, SNLS, HST and a compilation of nearby SN la) as described in [5]; b) 
observations of baryon acoustic oscillations in the SDSS Luminous Red Galaxy sample [8] , 
in the WiggleZ Survey [7] and 6dF Galaxy Survey [19] datasets; c) measurements of cos- 
mic microwave background temperature anisotropy from WMAP7 [20, 21]. In the following 
subsections, the statistical analysis used for those observables is described. 

3.1 Type la supernovae 

In Ref. [10], the g-parametrization given by eq. (2.1) had been investigated by using two 
different SN la datasets: the Goldl82 sample [22], processed with the light-curve fitter 
MLCS2k2 [25], and 115 SN la from the SNLS survey [23], analyzed with SALT [24]. In this 
work we make use of only one SN la dataset, the 288 SN la compilation from Kessler et al. 
[5] (sample "e" in their paper), and analyze it with both light-curve fitters MLCS2k2 [25] 
and SALT2 [26]. In fact, with SALT2 we use only 282 SN la, excluding the supernovae d086, 
90O, 93B, 3256, 930 and Gilgamesh since they gave negative values for a"^ (see eq. (3.17) 
for definition). 

A key quantity in SN la investigation is the distance modulus 

^JLth{z■, h, 9) := 5 log [Dl{z; 9)] + fioih). (3.1) 

Here h is the Hubble constant in units of 100 km/s Mpc^^, 9 denotes the set of cosmological 
parameters of interest other than h, 

;Uo(M := 5 log j' 10'^/0^"Vs) \ ^ 42.38 - 5 log h. (3.2) 

and 

«.(.;/., 9) = (l + ^)fj^i;4^, (3.3) 

is the dimensionless luminosity distance (in units of the Hubble distance today). In the 
following we discuss how to obtain constraints on cosmological parameters by using MLCS2k2 
and SALT2 output as our data. 

3.1.1 MLCS2k2 

In the MLCS2k2 fitting model [25], the variation among SN la light curves is described with 
a single parameter (A). Excess color variations relative to the one-parameter model are 
assumed to be the result of extinction by dust in the host galaxy and in the Milky Way. The 
MLCS2k2 model magnitude, observed in an arbitrary filter Y, is given by [5] 

where e is an epoch index that runs over the observations, Y' G {U,B,V,R,I} is one of 
the supernova rest-frame filters for which the model is defined, A is the MLCS2k2 shape- 
luminosity parameter that accounts for the correlation between peak luminosity and the 
shape/duration of the light curve. Furthermore, the model for the host-galaxy extinction 



is XyT^ = Cy',€{(1'Y' + bY'/Rv)Av, where Cy,e := XyTl/XyTQ, and ay', hy' are constants; 
as usual, Ay is the V band extinction, at B band peak (ay = 1, hy = 0), and Ry '■= 
Ay/E{B — V), the ratio of V band extinction to color excess, at B band peak. Finally, Xy ^ 
is the Milky Way extinction, Ky/y,, is the i^-correction between rest-frame and observer- 
frame filters, and ^ is the distance modulus. The coefficients My^^, py^e, and qy^e are 
model vectors that have been evaluated using nearly 100 well observed low-redshift SN la as 
a training set. Above, e = corresponds to the B band peak magnitude epoch. 

Fitting the model to each SN la magnitudes, usually fixing Ry gives /x. A, Ay and Iq, 
the -B-band peak magnitude time. 

In this work, to obtain the cosmological parameters constraints, we adopt the likelihood 
approach as described in [27]. Assuming that all SN la events are independent, for the 
MLCS2k2 data, we can write the likelihood as 



^ 1 
L{h,e,a'-') = ll ^exp 



'2^MLCS2k2(^'^i^ ) 
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where 
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In the expression above N is the total number of supernovae (A'^ = 288 for the MLCS2k2 anal- 
ysis), Hj is the j-th SN la distance modulus (u^ is its statistical uncertainty, both estimated 
from the MLCS2k2 fitter), Zj is the spectroscopically determined SN la redshift and 

Following [5], we assume an empty universe to project the redshift uncertainties onto distance 
modulus such that 

We checked that using an expression taking into account the cosmology has negligible effect 
on the constraints. Here, a^^ is the j-th SN la total redshift uncertainty 

(a^.)' = (^.%c)' + <„ (3.9) 

where a spec is the error due the redshift determination and apec = 0.0012 is an additional 
error due to peculiar velocities. In eq. (3.7), o"*"* is the intrinsic dispersion which is added 
in quadrature to the distance modulus and redshift dispersions. It represents intrinsic SN la 
dispersion (that are expected to exist even after corrections) and any potential systematic 
errors that have not been accounted for. In this work it is considered a free parameter to 
be determined from the data together with the other parameters ([27], see also [28]). The 
approach we follow is different from the usual MLCS2k2 x^ oiig in which o"*"* is estimated 
through an iterative procedure in which the x^ P^r degree of freedom is made unity for the 
Hubble diagram constructed from the nearby SN la sample [5, 29]. 

3.1.2 SALT2 

The aim of the SALT2 method [26] is to model the mean evolution of the SN la spectral 
energy distribution (SED) and its variations with a few dominant components. It makes use 



of a two-dimensional surface in time and wavelength that describes the temporal evolution 
of the rest-frame SED. When using the SALT2 light-curve fitter, the SN la rest-frame flux 
at wavelength A is modeled by the following functional form 

-^{t, A) = xo [Mo{t, A) + xiMi(t, A)] exp[c C(A)] (3.10) 

where t is the rest-frame time since the date of maximum luminosity in B-band and A is 
the wavelength in the rest-frame of the SN la. Mo(t, A) is the average spectral sequence 
whereas Mi{t,X) is an additional component that describe the main variability of SN la 
(more additional components might be added). C(A) represents the average color correction 
law [26]. In contrast to Mq, Mi and C that are properties of the global model (derived using 
a training set of SN la [26]), xq, xi and c are parameters associated with each individual 
supernova. 

To compare with photometric SN la data, the observer frame flux in passband Y is 
calculated as 
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where T (A) defines the transmission curve of the observer- frame filter Y, and z is the 
redshift. Each SN la light curve is fitted separately using eqs. (3.10) and (3.11) to determine 
the parameters xq, xi, and c with corresponding errors. The parameter c is a measure of the 
SN la color, xi is related to the stretch of the light curve and the peak rest-frame magnitude 
in the B band is given by 

m*B = -2.51og[xo f dX'Moit = 0, A')T^(A')]. (3.12) 

The j-th SN la distance modulus is modeled by 

Hj = m*^- — M + axij — f3cj , (3.13) 

and its statistical uncertainty by 

(aff = [af^f + a\aff + p\a^f + 2a(a7^"^) - 2/3(^7^') - 2aP{af') . (3.14) 

Here the quantities {a- ^)^, (o"|^)^, (crp^, a-^ \ a- ^ and a^-^^ are the components of the 
covariance matrix of (m^, xi, c). We remark that ttt.^, xi and c are derived from the fit to the 
light curves while the global parameters a, /3 and the absolute magnitude M are estimated 
simultaneously with the cosmological parameters and are marginalized over when obtaining 
the cosmological parameter constraints. 

Again assuming that all SN la events are independent, the likelihood for SALT2 data 
is given by 
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L{e, a, /3, 7W(/i, M), a^"*) = TT ^-= exp 
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where 



xIalt2{0, a, t3,M{h, M)y-') = J^ '^' ^^^^^^^'"''^J , (3.16) 



.^1 ''T, 



M{h,M) := fio{h) + M, Hj is given by eq. (3.13), 

4, = K)' + (<)' + (^^"*)' (3.17) 

and, for SALT2, o"^ is given by eq. (3.14) and N = 282. Alternatively, we can define the 
apparent magnitude in B band as the theoretical estimator for SALT2, 



BtJ^i! ^' a,P,M)=M- axij + /3c, + Slog [Dl(2,; 6>)] (3.18) 
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where ^j = {xij, Cj, Zj) denotes the set of observables and, for the considered model, 6 = 
{qi,qf (or qo),zt,T). 

With the above definition we can rewrite (3.16) as 

J— 1 J 

For the SALT2 data analysis we do not assume a fixed value for o"*"* [29] nor obtain its 
value by imposing, through an iterative process, that x^ P^r degree of freedom equals unity. 
As in the MLCS2k2 case, o"*"* is considered a free parameter and is determined from the data 
together with the other parameters [27]. Although, in this work, our marginalization proce- 
dure is purely numerical, in the Appendix we discuss how to marginalize (3.15) analytically 
over the parameter M. and (3.5) over /xq. 

3.2 Baryon acoustic oscillations and cosmic micro^vave background 

We start by defining the comoving sound horizon at the photon-decoupling epoch 

, , c /■V(i+^.) da 
J. i-^ \ — / (2, 20) 

^/^JQ a2/f(a)^l + (3[7M/4[7^o)a' 

where VL^q and ri{,o are, respectively, the present value of the photon and baryon density 
parameter. In eq. (3.20) z* is the redshift of photon decoupling and is well approximated by 
the formula given in [30]. In accordance with WMAP7 [21], in this work we use z^^ = 1091 
exactly. We checked that variations within the uncertainties about this value do not alter 
significantly the results [29]. Another relevant quantity is the redshift of the drag epoch 
{zd ~ 1020), when the photon pressure is no longer able to avoid gravitational instability of 
the baryons. 

To obtain the BAO/CMB constraints we use the "acoustic scale" 

i. = >r^ , (3.21) 

where ^^(-z*) = c L* dz' /H{z') is the comoving angular-diameter distance. Percival et al. 

[8] measured rs{za)/Dv{z), at z = 0.2 and z = 0.35. Here, Dy{z) := [d\{z)cz/H{z)'^ is 
the "dilation scale" introduced in [31]. In [19] the 6dF Galaxy Survey also reported a new 
measurement of Vs/Dy at z = 0.106 and, more recently, the WiggleZ team [7] obtained results 
at z = 0.44, z = 0.60 and z = 0.73 (see table 1). Combining these results with the WMAP7- 
year [21] value I a = 302.44 ± 0.80 give the values shown in table 1 for r.'^f^''*^ . ^H. By 

J ^ ' ^ O Dv{ZbAO) '•s(2*) •' 

using the WMAP7 [21] recommended values for rs{Z(i) and rs{z^:) we obtain rs{zii)/rs{z^) = 
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32.35 ±1.45 


18.34 ±0.59 


10.56 ±0.35 


8.82 ±0.68 
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5.70 ±0.31 


Dvl^BAo) r,(z,) 
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30.95 ±1.46 


17.55 ±0.60 


10.11 ±0.37 


8.44 ± 0.67 


6.69 ±0.33 


5.45 ±0.31 
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Table 1. Values of „'';^^''^ , [7, 8, 191, „'^f<^-^ , ^^ and n'^f^^"^ , for different values of 
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1.045 lb 0.016. By inserting this ratio in the expression above we obtain the BAO/CMB 
constraints ( r, f. . n ) also exhibited in table 1. 
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Now we can write the x^ for the BAO/CMB analysis as 
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and 
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is the inverse covariance matrix derived by using the above results together with the cor- 
relation coefficients r = 0.337, r = 0.369 and r = 0.438 calculated for the Tg/Dy pair of 
measurements at z = (0.2,0.35), z = (0.44,0.6) and z = (0.6,0.73), respectively [7, 8]. 

4 Results 

Through all this work, to obtain the probability distributions (PDFs), the Metropolis-Hasting 
algorithm has been used [32]. Generally, to obtain the PDFs, 5000 chains were generated 
with 10000 points for each chain. We first consider models in which the final value of the 
deceleration parameter is fixed to g/ = —1. Cosmological models in this special class asymp- 
totically tend to a de Sitter phase in the future. Several popular cosmological models like 
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Figure 3. 68% and 95% confidence contours in the plane arctan(zt) vs. arctan(T), for qi — 0.5 and 
Qf = —1, for SN la data only, for both MLCS2k2 (in red) and SALT2 (in green). We marginalized 
over all other parameters with flat priors. The dots indicate the best-fit values. 



flat ACDM, quartessence Chaplygin, DGP etc, belong to this class as discussed in section 
2.2. 

By fixing qj = —1, we are left with only two free model parameters, the transition 
redshift zt and its duration r. In figure 3 we display 68% and 95% confidence level (C.L.) 
regions in the (arctanzijarctanr) plane allowed by supernovae experiments. As mentioned 
in section 3, in this work we are considering the Kessler et al. [5] SN la compilation. From 
figure 3 we can see that, at high confidence level, zj ^ is not permitted, indicating that, 
from SN la observations, a transition occurred in the past. Furthermore it is also clear from 
the figure that, SN la observations alone are not able to strongly constrain the parameters zt 
and T. For instance, a high value of the transition redshift (let say zt > 5) is allowed by SN 
la observations if the transition is slow (r > 1). Since SN la observations probe the universe 
only up to z ~ 1 — 2, the distance to an object, say, at z < 1 in a universe with r > 1 and 
high Zt, can be similar to the distance to the same object, in a universe with zt < 1 but with 
a faster transition (smaller r) and this explains the shape of the contours. 
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Figure 4. 68% and 95% confidence contours in tlie plane zt vs. r, for qi ~ 0.5 and qf — —1, 
for BAO/CMB data only. We marginalized over all other parameters with flat priors. Left: result 
obtained with only the two measurements from [8]. Right: result obtained with all six measurements 
considered in this work. The dots indicate the best-fit values. 



In figure 4 we show BAO/CMB constraints on zt and r in the case qf = —1. In the left 
panel we display 68% and 95% confidence regions obtained using only the two measurements 
of rs{zd)/Dv{z) obtained by Percival et al. [8]. In the right panel we display the results 
when six measurements of rs{zd)/Dv{z))[7, 8, 19] are taken into account. The improvement 
these new observations bring in constraining the model parameters is evident. Notice also 
that this test is complementary to the SN la one and by combining them we can considerably 
reduce the allowed parameter space. 

Since the SN la data are independent from BAO/CMB we can write the combined x^ 
statistics as 



X 



-2ln Lsn + Xbao/cmb 



(4.1) 



where x'bao/cmb ^^ given by eq. (3.22) and Lsn by eq. (3.5), when adopting MLCS2k2, or 
eq. (3.15) when the SALT2 light-curve fitter is used. In figure 5 we show the 68% and 95% 
C.L. regions in the {zt,T) plane imposed by the combined data sets. In the figure the dotted 
lines represent the ACDM (r = 1/3) and the DGP (r = 0.5) models. The figure indicates 
that the combination of the BAO/CMB data set with SN la analyzed with MLCS2k2 tend 
to favor the DGP model, while the combination of BAO/CMB plus SN la analyzed with 
SALT2, tends to favor the ACDM model. 

We now consider the situation in which we do not impose a fixed value for qj. In this 
quite general case, we would have to determine three model parameters, namely zt, r and 
qf (we still assume qi = 1/2). However, instead of working with qf, we find more convenient 
to use qo, the present value of the deceleration parameter. By using eq. (2.11), one can 
re-express q{z) in terms of qo (see eq. (2.12)). We remark that possible r values are not 
arbitrary but constrained by eq. (2.13). 

In figure 6 we display 2D marginalized (flat prior) C.L. regions (68% and 95%) in the 
{zt, r), {zt, qo) and {qo, zt) planes imposed by the combination of SN la (using MLCS2k2 and 
SALT2 fitters) and BAO/CMB tests. The marginalized (fiat prior) ID likelihoods (nornial- 
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Figure 5. 68% and 95% confidence contours in the plane zt vs. r, for qi = 0.5 and q^ = —1, for 
SN la (MLCS2k2 in red and SALT2 in green) combined with BAO/CMB (six measurements). We 
marginahzed over all other parameters with flat priors. The dashed lines denote the values of the 
parameters that correspond to the models ACDM (r = 1/3) and DGP (t = 1/2). The dots indicate 
the best-fit values. 





MLCS2k2+BAO/CMB 


SALT2+BA0/CMB 


Zt 


^ .^+0.13 (+0.45) 
'-'•^"-0.10 (-0.18) 


n fizL+O-13 (+0-30) 
^•"^-0.07 (-0.12) 


T 


+0.16 (+0.27) 
'-'•^'-0.20 (-0.40) 


p, oc+0.11 (+0.19) 
^•"^"-0.17 (-0.30) 


qo 


p, oi+O.ll (+0.20) 
^•"^-^-0.11 (-0.23) 


P, cro+0.17 (+0.28) 
^•^"^-0.13 (-0.27) 


O'int 


^ , o+O.Ol (+0.02) 
^•-"^■^-O.Ol (-0.02) 


p, ,„+0.01(+0.02) 
^•-^^-0.01 (-0.02) 


a 


- 


n 004+°°^^ (+0.022) 
u.u»'i_g Q-^-^ (_0.022) 


/3 


- 


„ ,j,+0.10(+0.21) 
'^•-^'^-0.11 (-0.21) 



Table 2. Summary of the best-fit values for all parameters when using SN la + BAO/CMB, including 
the 68% and 95% confidence intervals. For each parameter we marginalized over all other parameters 
with flat priors. 



13 




o 'f 

B^CMB 



Figure 6. 68% and 95% confidence contours in tlie planes zt vs. t (bottom left), zt vs. go (middle 
left) and go vs. r (bottom center), considering qi = 0.5, for SN la (MLCS2k2 in red and SALT2 
in green) combined with BAO/CMB (six measurements). The dashed lines denote the values of the 
parameters that correspond to the models ACDM (t = 1/3) and DGP (r = 1/2). We also show the 
ID PDFs [C/ Cmax) for Zt (top), q^ (middle center) and r (bottom right). In all cases we marginalized 
over all other parameters with flat priors. 

ized to its maximum values - C/ Cmax) for each cosmological parameter, zt, t and qo, are also 
shown. In the figure are also displayed dashed lines corresponding to the flat DGP (r = 1/2) 
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Figure 7. 68% and 95% confidence surfaces in the space (go, ''', zt), for qi — 0.5, for SN la (MLCS2k2 
in red and SALT2 in green) combined with BAO/CMB (six measurements). We marginahzed over 
all other parameters with flat priors. The lines represent the flat models ACDM and DGP. 





08 0.09 0.10 0.11 0.12 0.13 0.14 0.15 0.16 



Figure 8. ID probability distributions (C/Cmax) for o'mt, marginalized over all other parameters 
with flat priors. The shaded regions indicate the 68% and 95% confidence intervals. 



and ACDM (r = 1/3) cosmological models. In table 2 we show the best-fit and the do7o 
and 95% limits on several quantities after marginalizing over all the other parameters with a 
flat prior. In figure 7 the 68% and 95% C.L. 3D regions for the parameters qo, zt and r are 
shown. The curves corresponding to the flat ACDM and DGP models are also displayed. It 
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Figure 9. Bottom left: 68% and 95% confidence contours in tlie plane a vs. /3, considering SALT2 
+ BAO/CMB and marginalizing over all other parameters wilit flat priors. Top: ID probability 
distribution {C/Cmax) for a. Bottom right: ID probability distribution {C/Cmax) for /?• 



is clear from the figure that the flat ACDM model is favored when we combine BAO/CMB 
with SN la analyzed with SALT2, while flat DGP is preferred if we use the MLCS2k2 fitter. 
In figure 8 we show the ID probability distribution, C/Cmax (marginalized with a flat prior), 
for the the intrinsic dispersion o"j„j and in figure 9 the 68% and 95% C.L. regions in the (a, (3) 
plane are displayed together with ID distribution (C/Cmax) for these parameters. 



5 Discussion 

In this work we considered recent SN la, BAO and CMB observations to investigate the 
transition from cosmic deceleration to acceleration. Instead of adopting the traditional x^ 
treatment, the SN la analysis was performed with the likelihood method, as described in [27]. 
Following [10] we considered a kinematical approach using a kink-like parametrization for the 
deceleration parameter. We first considered the simple case in which the models have a final 
de Sitter phase such that qf = —1. Since [10] focus on this case, it is worthwhile to compare 
the results obtained in the present work with those in that paper. First let us remind that 
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in [10] two different SN la datasets were used: tfie Goldl82 sample [22], processed with tlie 
ligfit-curve fitter MLCS2k2, and the SNLS survey [23], analyzed with SALT. Here we make 
use of only one SN la dataset, the Kessler et al. [5] (sample combination "e") compilation, 
and analyze it with two light-curve fitters (MLCS2k2 and SALT2). In [10], the BAO/CMB 
constrains upon the model parameters were obtained using only two BAO measurements, 
whereas here we have used six. There is no difference between what has been referred to in 
[10] by the Sk/Dy test, and what is denominated here the BAO/CMB cIa/Dv test. Our 
figure 4 (left panel) is essentially the same as figure 1 (right panel) of [10]. It is worth 
mentioning that dA{z^)/ Dv{zbao) is equal to zbao times the ratio TZ{z^) / A{zbao) , where 
TZ is the CMB shift parameter and A is the BAO scale [31]. Notice that this ratio eliminates 
the dependence on ilmo present in the definition of these quantities and, therefore, is less 
model dependent. 

In [10] a tension was found between the results obtained using the Goldl82 dataset 
and those obtained with the SNLS dataset. The comparison was performed only for the 
special case qf = —1 and the origin of the tension was attributed to the datasets themselves. 
Here, with a more general approach (we do not fix qj), we showed that the origin of the 
tension might not be the dataset (we used only one) but the light-curve fitter adopted. 
For flat space, when combined with BAO/CMB, SALT2 tends to favor standard ACDM 
while MLCS2k2 favors the DGP model (see figure 7). This confirms, with a quite model 
independent approach, a similar result obtained in [29] by using an information-criteria model 
selection. We remark that the discrepancy between the MLCS2k2 and SALT2 light-curve 
fitters has originally been pointed out in [5]. Our results are also consistent with those of [5] 
that considered the more conventional ACDM and fiat tuCDM models. For instance, when 
combining MLCS2k2 with BAO and CMB measurements they obtained w = —0.76 as best- 
fit (SN la sample combination "e") for flat ifCDM, while for SALT2 they got w = —0.96, 
that is in good agreement with ACDM. Interestingly they found that this difference of ~ 0.2 
in the dark energy equation of state parameter is reduced to < 0.04 if only the SDSS or 
Nearby+SDSS samples are used. Including ESSENCE gives also consistent results. By 
considering the mean difference in distance modulus, A/i = ij-salt2 — IJ'MLCS2k2, as a function 
of redshift, for different SNIa samples, they showed that, unlike the other three samples, the 
mean A/u is considerably different from zero (and positive) for SNLS and HST samples. 
Therefore, including these two samples push w to smaller values, as compared to MLCS2k2, 
when adopting SALT2. Similarly, the same happens in our case. For instance, assuming 
qi = 1/2, qf = —1 and a fixed arbitrary value zt, it can be shown that distances are smaller 
for DGP (r = 1/2) than for ACDM (r = 1/3) and this qualitatively explains our results. We 
finally remark that a comprehensive discussion about possible sources for the discrepancy 
between MLCS2k2 and SALT2 is performed in [5] that trace it to rest-frame ultraviolet 
modeling and also to the intrinsic color and dust extinction treatment in both fitters. 

The main results of this work are displayed in figures 6, 7 and table 2. It is now clear that 
progress in determining the transition redshift and/or the present value of the deceleration 
parameter, depends crucially on solving the issue of the difference obtained when using these 
two light-curve fitters. This points out to unresolved systematics as emphasized in [5, 29]. 
Indeed, in comparison, answering other questions like which parametrization is more adequate 
to determine qo or zt, is currently less important. 

A relevant issue, not considered in this work, is to quantify how universal is the as- 
sumption qi = 1/2. We are currently investigating how good this hypothesis is by taking 
into account more observables. Another interesting point that deserves further study is how 
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well future SN la, BAO and CMB surveys will able to constrain the model parameters. The 
results of these investigations will be published elsewhere. 
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A Analytic Marginalization over /jq in MLCS2k2 and over A^ in SALT2 
SN data analysis 

A.l Considering MLCS2k2 output 
A. 1.1 Uncorrelated observations 

Writing explicitly the dependence of x'MLCS2k2 ^^ f^o ™ ^^- ("^-6)) we have 

XMLCS2k2(/^0, ^, CTint) = }, ^ 2~? — ''~\ ' (^•^) 

^ crrp^iicTint) 

Introducing the following auxiliary quantities 

Bie,a.^,):=Y^ ^^^-'':f^['^-^'^^\ (A.3) 

^ 1 

C{a,nt):=^^-, r, (A.4) 

^ 1 

Dia,nt) ■■= n } w^r ' (A-5) 

we can rewrite eq. (A.l) as 

XMLCS2k2(M0, 0, CTint) = ^{6 , Uint) - 23(6, aint)fJ-0 + C'(c7i„t)/io, (A.6) 

consequently, the likelihood, eq. (3.5), will be expressed by 



L(//o, 0, (Tint) = D{aint) exp 
The marginalized likelihood will be 



^O^CFint) . TjlQ \ C{(Tint) 2 
+ B{e, aint)fJ^O 7, Mo 



/ 2tt 
L*{6,aint) = D{aint)\ yrr r exp 
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2C{a, 



int) 



(A.7) 



/ + 00 
L{fio,e,aint)dfio. (A.8) 

-oo 

Substituting eq. (A.7) in (A.8) we obtain 

' Aie,a^nt) , B^ie^aintY 



(A.9) 



Finally, we can write 



*2 



X*MLCS2k2{^J^int) ■— "2 In L*(0, <Tmi) — A{6 , aint) 



A. 1.2 Correlated observations 

In the case of correlated observations we have 

1 



B^e,a, 



int) 



C{(Tint) 



2 In 



D{aint)] 



27r 



C{aint 

(A.IO) 



L(/xo,0,cr, 



int) 



exp ( --X^C-'X ] , 



V(27r)^detC 
where, for MLCS2k2 data, 

Xi = [fii - fith{zi;fio,0)] = [/ij -5log'DL{zi;e) - hq] 
Following [4] we will write the off-diagonal terms of C as 



fc=i 



(A.ll) 
(A12) 

(A13) 



where the sum is over M systematics Sk and AS^ is the uncertainty associated to Sk- Given 
this expression, we find that C* do not depend on /ig- 



bitroducing a new vector 



such that 



we can write 



Xi = [m - hlogV Lizi'.e)], 
X = X -/iol, 



x' = x^c-^x = {x- ^i^ifc-\x - /.oi). 

Introducing new auxiliary quantities 

a{0,(Tint) 

b{e,aint) 
e{aint) 
d{aint) 



x^c-^x, 

x^c-^i, 



we have 
and 



V(27r)^detC' 

XMLCS2k2(/^0> ^5 <7int) = a{6, aint) — 26(0, (Jint)lJ'0 + e{aint)fJ', 



L(/xo, 0, aint) = d{aint) exp 



a{0,aint) , ,,„ ^ e{aint) 2 
+ b{0,aint)fio — /^o 



Finally, marginalizing over fiQ, we have 



and 



L* {6, aint) — d{aint)\ 



X* MLCS2k2{^^(^int) - a{6,aint) 



27r 



— -exp 

e[aint) 

b^ {9, aint) 
e{aint) 



a{9,aint) b'^{e,aint) 



'^e{aint) 



2 In 



d{ai, 



nt 



I 2tt 
ei<Jint) 



(A.14) 
(A15) 
(A.16) 

(A.17) 
(A18) 
(A19) 
(A.20) 

(A.21) 

(A.22) 

(A.23) 
(A.24) 
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A. 2 Considering SALT2 output 
A. 2.1 Uncorrelated observations 

Let us write eq. (3.19) as 



XsALT2^>a,P,-M,<7mi := > 3-7 o ^ : 



where 



2 f a \ 2 |/u\2|22 io2 2|2 

Introducing the following auxiliary quantities 



B [6, a, /3, aint) ■= 2^ 



0"'p j(<7mi) 



N 



C'{a,/3,aint) ■■= ^ 



1 



TV 



D' {a, (3, aint) := J] 



a^^{a,P,aint)' 



1 



(A.25) 
(A.26) 

(A.27) 
(A.28) 
(A.29) 
(A.30) 



(TT,i(a,/3,o-jnt)V27r' 
we can rewrite eq. (A.25) as 

XsALT2(^, a, (3, M, aint) = A' {6, a, /3, aint) - 2B'{0, a, (3, a,nt)M + C'{a, f3aint)M^ (A.31) 
consequently, the likelihood will be expressed by 



L{6,a,P,M,aint) = D'{a, f3,aint) exp 
The marginalized likelihood will be 



A' {6, a, ^, aint) , o/^a „. a _ na^ C"(a,/3, cTmi) ^ v2 



B'{e,a,P,aint)M 



/+00 
L(6>,a,/3,A^,cJi„t)o!A^. 
-00 

Using eq. (A. 32) in (A. 33) we have 



2 

(A.32) 



(A.33) 



-M' 



2tt 



L* {6, a, /3, aint) = D'{a,/3,aint)\ -^^ — ^ r exp 

y C'{a,f3,aint) 



Finally, we can write 



h<2 



A' {6, a, p, aint) B'\e,a,P,aint) 
2 2C'{a,^,aint) 



A' {e, a, 13, a, 



intj 



X*sALT2(^)a>/3,o-mt) := -21nL*(0,a,/3,crj„t) 
2 In D' {a, (3, aint)] 



(A.34) 
(A.35) 



B''{0,a,(3,aint) 



C'{a,/3,aint) 



27r 



C'{a,l3,aint) 
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A. 2. 2 Correlated observations 

In the case of correlated observations we just need to choose 

Xi = [m*s^i-mth{0,a,^,M,Zi)] = [m*^^^ - 5logVL{e, Zi) + axi^t - (3ci - M]. (A.36) 
Following [4] again, we will write the off-diagonal terms of C as 



M 



a 



jLj _ sr^ dmth{0,a,f3,M,Zi) dmth{9,a,l3,M,Zj) 



u 



E 

fe=i 



dSk 



dSk 



(ASk 



(A.37) 



Now we find that C*- do not depend on ^A. 



V 



Introducing a new vector 

Xi = [m*^^ - 5log'DL{6,Zi) + axi^i - (3ci], 

such that 

X = X -Ml, 

and introducing the auxiliary quantities 

a'{e,a,f3,aint) 
b\e,a,f3,aint) 

e'{a,/3,aint) 

d'{a,(3,ai„t) : 



x^c-^x, 

x'^c-H, 

fc-% 

1 



(A.38) 
(A.39) 

(A.40) 

(A.41) 
(A.42) 
(A.43) 



v/(2^)^detC 
we have 

XsALT2(^> a, l3,M,aint) = a {6, a, /?, aint) - 2b' {0, a, f3, aint)M + e'(a, /3, ai„t)X^ (A.44) 
and 



L{e,a,P,M,aint) = d'{a, /3,aint)exp 
Finally, marginalizing over Ai, we have 
L*{e,a,l3,aint) = d'{a, /3,aint'' 



a'{6,a,/3,aint) , ,, 



+ b'{e,a,f3,aint)M 



e {(Tint, 



M' 



(A.45) 



27r 



e'{a,j3,aint, 



exp 



a'{6,a,P,aint) b''^{0,a, (3,aint) 
2 2e'{a,P,aint) 



(A.46) 



and 



X*5ALT2(^> "; P, CTint) = a' {8, a, /?, a, 



int) 



b"{e,a,/3,aM) 
e'{a,/3,aint) 



-2 In 



d'{a,(3,aint)\ - 



2tt 



d {a, 13, a int 

(A.47) 
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